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ABSTRACT 

Ultraluminous infrared galaxies (ULIRGs) have several types according to dominance of starburst or AGN 
component. We made stellar population analysis for a sample of 160 ULIRGs to study the evolution of 
ULIRGs. We found that the dominance of intermediate-age and old stellar populations increases along the 
sequence of Hll-like ULIRGs, Seyfert-HII composite ULIRGs, and Seyfert 2 ULIRGs. Consequently the typi- 
cal mean stellar age and the stellar mass increase along the sequence. Comparing the gas mass estimated from 
the CO measurements with the stellar mass estimated from the optical spectra, we found that gas fraction is 
anti-correlated with the stellar mass. Hll-like ULIRGs with small stellar masses do not possess enough gas 
and the total mass, and therefore have no evolution connections with massive Seyfert 2 ULIRGs. Only massive 
ULIRGs may follow the evolution sequence toward AGNs, and massive Hll-like ULIRGs are probably in an 
earlier stage of the sequence. 

Subject headings: galaxies: infrared — galaxies: evolution — galaxies: starburst — galaxies: stellar content 



1. INTRODUCTION 

Ultraluminous infrared galaxies (ULIRGs) are advanced 
mergers of gas-rich galaxies. They have an infrared lu- 
minosity greater than lO^'^ Lg, in 8-1000 pL m band (e.g. 
ISanders & MirabeilfT99l iLonsdale et al.\Mm . Their high 
infrared luminosity comes from violent starbursts with signif- 
icant contributi on from active ga lactic nuclei (AGNs) in some 
ULIRGs (e.g. F arrah et al.ll2003h . 

ULIRGs can be classified according to dominance of star- 
burst or AGN component, namely as the Hll-like ULIRGs, 
LINER ULIRGs, Seyfert-HII com posite ULIRG s Type 1 
ULIRGs and Seyfert 2 ULIRGs. ISanders et all (Il988allbh 
proposed an evolutionary scenario of ULIRGs from cool- 
ULIRGs to warm-ULIRGs and finally to QSOs. The cool- 
ULIRGs have f 25 /.fen < 0.2 an d are preferably star-forming 
galaxies (iHeckman et alJ [T987 t). i.e. the Hll-like galaxies. 
Here /25 and /eo are flux densities at 25 fim and 60 /im, 
respectively. The warm-ULI RGs have /25//60 > 0-2 and 
prefer ably host an AGN (de Grijp et al.ll985tlBeichman et all 
Il986l) . ULIRGs are important to understand the evolu- 
tion of merger galaxies and the possible connections be- 
tween the starburst an d AGNs. Fu rther i maging and spec- 
troscopy studies (e.g. Surace et al. 1998; Surace & Sanders 
[1999; Surac e et al. 2000; Lfpari et al. 2003; Schw eitzer et al.. 
I2PO6; Netz er et a 



I2OOI; (\^iUeux et al.l 2009) and numerical 



simulations (e.g. iKor mendv & Sanderill992; .Springel et al.l 
UpOS; Naa b et alj 1200 6) support this evolutionary scenario. 
Ho wever, this s tandar d evolution sequence is que stioned 
by IColina et"an (1200 Ih: [Genzel et al (2001 ); Tacc oni et alJ 
(I2002h : [Rodriguez Zaurin et al.l (l2010l) . 

If the standard evolutionary scenario is true, ULIRGs 
should experience the continuing conversion of gas compo- 
nent to stellar populations in their lifetime. ULIRGs in the 
earlier evolutionary stage are expected to possess more gas 
and young stars, less old stellar populations, smaller stellar 
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mass and larger gas mass fraction, compared with those of 
ULIRGs in a later evolutionary stage. There should be enough 
gas in early stage ULIRGs to be converted into stars when 
they evolve to a later stage for a larger stellar mass. 

The stellar population analysis of ULIRGs can shed light on 
ULIRG evolution. Previous s uch analys es were only carried 
out on very smafl samples (e.g.lCanalizo & Stockton 2000b A 
20011: [Rodrigu ez Zaurin et al.l l2007l l2008i; iSoto & Martini 



2010h . A slightl y larger sample of 36 local ULIRGs was stud- 
ied by Rodrigue z Zaurin et al.l (12009. 2010). In contrast to the 
standard evolutionary scenario, their stellar population analy- 
sis based on the optical spectra exhibits no significant differ- 
ences of the stellar ages among the Hll-like ULIRGs, LINER, 
and Seyfert 2 ULIRGs. Note th at the physical para meters of 
ULIRGs ai-e rather scatter (e.g. 'Veifleux et al.l 120091) . there- 
fore, the possible evolutionary of ULIRGs can not be revealed 
if only a small sample is applied to. 

In this work, we analyse the spectra of 160 ULIRGs 
available from the Sloan Digital Sky Survey (SDSS, 
lAbazajian et al.l2 009^ for stellar populations. In Section 2, we 
introduce our ULIRG sample. The stellar population analysis 
of the sample is presented in Section 3. In Section 4, we dis- 
cuss the inferred stellar age, mass and gas for different types 
of ULIRGs. Several factors which affect our analysis are dis- 
cussed in Section 5, and conclusions are given in Section 6. 

In this paper, we adopt Hq = 70 km s^^ Mpc^^, flm = 0.3 
and JIa = 0.7. 

2. THE ULIRG SAMPLE 

For reliable stellar population analysis, we work on a large 
sample of ULIRGs which have optical spectra covering the 
wavelength regions of 4000 A break and strong stellar ab- 
sorptions near 4000 A. 

The SDSS database (DR7. [Abazaiian et alj[200^ released 
the spectra of more than 920,000 galaxies in the range 
of 3800 A to 9200 A in the observer frame . Follow- 
ing the procedures described in iHou et akl (l2009l) . we found 
398 ULIRG candidates from the SDSS DR7 specti-al cata- 
log (Abazaiian et al. 2009) and IRAS Faint Sources Catalog 
( iMoshir et al. 1992). We checked them from NASA/IPAC Ex- 
tragalactic Database (NED) and SDSS images. We found that 
38 of them are wrong identifications of ULIRGs due to either 
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Fig. 1 . — BPT diagrams for line ratios of 185 narrow-line ULIRGs. 
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Fig. 2. — Redshift distribution of 160 narrow line ULIRGs, after we dis- 
carded 22 ambiguous objects and 3 Seyfert-HII composite ULIRGs with bad 
pixels from 185 narrow-line ULIRGs. 



the known redshift of a ULIRG in the NED not consistent with 
the redshift from the SDSS database or an obviously nearby 
large bright galaxy as an IRAS source but with a mis-matched 
SDSS redshift. The spe ctra of 360 identified ULIRGs are 
then fitted and classified (iHou et al.ll2009i) . Among them, 73 
ULIRGs have broad lines and are of Type 1, and 287 ob- 
jects are narrow-line ULIRGs. We will not work on Type 1 
ULIRGs for stellar population studies in this paper because 
their spectra are dominated by strong AGN features (e.g. the 
power-law continuum, the strong and broad emission lines, 
the strong Fen features). 

Based on the line ratios in the BPT diagrams (see FigurelT] 
Haldwin et al. 1981; Kewlev et al. 2006), 185 of 287 narrow- 
fine ULIRGs are further classified as 32 Hll-hke ULIRGs, 
77 Seyfert-HII composite ULIRGs, 48 Seyfert 2 ULIRGs, 6 
LINER ULIRGs, and 22 ambiguous objects. We omitted the 
other 102 narrow-line ULIRGs from stellar population anal- 
ysis, because one or more lines either are redshifted outside 
the wavelength range or do not have a good signal-to-noise 
ratio for the classification in the BPT diagrams. We also dis- 
carded 22 ambiguous ULIRGs for our analysis because they 
are shown as one type in one BPT diagram but another type in 
the other diagram(s). Three Seyfert-HII composite ULIRGs 
are further skipped because more than half of pixels in their 
SDSS spectra are masked as bad. Therefore, we will study 
the stellar populations of 160 well-classified ULIRGs. Their 
redshift distribution is shown in Figure |2] 

In Figure |3] we show the combined spectra for four types of 
narrow-line ULIRGs. The spectrum of every ULIRG is first 
normalized by the median flux in the wavelength 4010 A - 
4060 A. The combined spectra are the weighted mean of the 
normalized spectra according to their signal-to-noise ratio in 
the wavelength range of 5300 A -5800 A in the rest frame. 
The combined spectra clearly show that along the sequence 
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Fig. 3. — Combined spectra of four types of narrow-line ULIRGs. 



of Hll-hke, Seyfert-HII composite, Seyfert 2, and LINER 
ULIRGs, the normalized continuum fluxes below 4000 A de- 
crease and the fluxes above 4000 A increase, which roughly 
indicates that the fraction of old stellar populations increases 
along the sequence. 

3. STELLAR POPULATION ANALYSIS 

The stellar population analysis is an important tool to 
study star formation in galaxies and evolution of galax- 
ies. It has been applied to various type of galaxies, 
e.g. HII galaxies (e.g. Schmitt et al. 1996; W estera et alJ 
12004 . AGNs (e.g ICid Fernandes et al.. ■2004a.bl) . infrared 



selected galaxies dChen et alJ l2009l l2010h. U LIRGs (e.g. 

m 



[Rodriguez Zaurin et al.il2010t iMeng et al.l l201(il) . Here, we 
use it for a large ULIRG sample to study the evolution se- 
quence of ULIRGs. 

We noticed that the fiber diameter of SDSS spectrograph is 
3", corresponding to 3 kpc for a galaxy at redshift z ~ 0.05 
or ~ 18 kpc at z ^ 0.25. For any galaxies with different red- 
shifts, the observed spectrum comes from central regions of 
diff erent size. We use the software, STARLIGHT ( version 
04. ICid Fernandes et al.l2005l:lMateus et al.l2006l;lAsari et al.l 
I2OO7I) . to analyze the stellar populations of ULIRGs in the 
central regions. We will discuss the aperture effect later 

The STARLIGHT fits an observed galaxy spectrum 0(A) 
with a combination of N simple stellar populations (SSPs). 

■* http://www.starlight.ufsc.br 
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Fig. 4. — One example for spectrum fitting. We show the normalized observed spectrum (i^A/^A4020' black) and the modeled spectrum (green) in the 
left-upper panel and the residual in the left-lower panel. The key parameters from the fitting: x^, mean stellar age, stellar mass and extinction parameter are 
given in the right side. The fractional contribution to the model flux (Xj) and the stellar mass (pj) of the j'th SSP are plotted in the right. 



The intrinsic extinction due to the foreground dust in the host 
galaxy is considered and parameterized by the V-band ex- 
tinction Av- The line-of-sight stellar motions are modeled by 
a Gaussian distribution G{vq, ct*) centered at the velocity vq 
and with the dispersion ct*. The modeled spectrum A/ (A) is 
described by: 

M(A) = A/aoE2;,6,(A) G{vo,<J,)]r{X) (1) 

where M\g is the synthetic flux normalized at Aq, Xj is 
the fractional contribution of the jth SSP to the model flux 
at Ao, bj{X) is the normalized spectrum of the jth SSP, 
and r(A) = 10~''-4[^(-^)^'^^ol is the reddening term. The 
best fitting to search for the minimum = J2xi(^W ^ 
A/ (A)) • w(A)]^, where w{X) is the weight factor, and w{\)^^ 
is the uncertainty of observed spectrum 0(A) given in the 
SDSS database. We used a base of TV 45 SSPs with 
different a ges and metallicities from the evolution synthesis 
models of iBruzual & Chariot! (I2003h. and adopted the initial 
mass function from IChabrier' ('2003') and t he Padova 1994 
evolutionary tracks (Alongi et al. 1993; Bre ssan et al.lll993t 
iFagotto et al. 1994; Girardi et al .1 11996) and STELIB library 
( |Le Borgne et al . 2003|). The base comprises of star popula- 
tion with 15 different ages between 1 Myr to 13 Gyr at each 
of the three metallicities: 0.2, 1, and 2.5 Z©. Here Z© is the 
metall icities of the Sun. The reddening law of ICalzetti et al.l 
(11994) is adopted in the fitting. 

We obtained the spectra with uncertainties of 160 ULIRGs 
from SDSS. We corrected the Galactic reddening effect, and 
converted each spectrum to the rest-frame. The STARLIGHT 
was used to fit the continuum and absorption features. The 
emission lines and sky lines are discarded in the mask file 
of STARLIGHT. The spectrum pixels without error measure- 
ments or with negative flux values are excluded. In addi- 



TABLE 1 

Mean fractional contributions of different stellar 
populations and the power-law component to the model 
flux at 4020 a for different types of ulirgs 



Type 


young 


intermediate-age 


old 


power-law 




« 10* yr) 


(10** - 10^ yr) 


(> 109 yr) 




Hll-like 


0.622 


0.261 


0.117 




Composite 


0.355 


0.408 


0.138 


0.099 


Syefert 2 


0.177 


0.500 


0.198 


0.125 


LINER 


0.112 


0.579 


0.261 


0.048 



tion, the Na D doublet, 5870 A - 5905 A, and other three 
bands, 6845 A - 6945 A, 7550 A - 7725 A, and 7165 A - 
7210 A are also mas ked, due to the bugs in SSP model 
(iLe Bo rgne et al."2003) or the large fitting residual in some re- 
gions iMateus et al. 2006). Considering the redshift coverage 
of our sample, we restricted the fitting in the wavelength range 
3400 A — 6700 A for each spectrum. A power-law component, 
F(A) cx A", is used in the STARLIGHT to account the AGN 
contribution to the observed continuum (Cid Femandes et alj 
2004b, 2005), and a = - 1.5 (Richstone & Schmidt 198^ 
Cid Femandes et al]|2004bl) is adopted. The random Markov 
Chains method was used in the STARLIGHT, which needs an 
input of integer see d. Different s eeds slightly affect the fitting 
results. Following Meng et al.l ( 1201 0). we fit each spectrum 
100 times with different seeds. The final fitted parameters are 
the mean values from 100 runs. One example of the spectrum 
fitting is shown in Figure |4] 

The stellar populations in a galaxy can be described by the 
fractions of young stars of tj < 10^ yr, intermediate-age stars 
of 10^ yr < t, < 10^ y r, and old stars of tj > 10^ yr 
dCid Femandes et al.ll2005l) . Here tj is the age of the jth SSP. 
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Fig. 5. — Fraction distributions of stellar populations for different types of 
ULIRGs from top to bottom. 
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Fig. 6. — Light-weighted mean stellar ages of ULIRGs estimated by the 
stellar population analysis are compared with the D„(4000) obtained by the 
MPA/JHU group. The solid line is the best fit for all ULIRGs. The distribu- 
tions of the mean stellar ages for each type of ULIRGs are shown in the right 
panels. 



We obtain the fractional contributions of these three stellar 
populations of different ages and the power-law contribution 
to the model spectrum flux at 4020 A. For each type of 
ULIRGs, we got their means as given in Table [T] Because 
there are only a small number (only 6) of LINER ULIRGs, 
we will not discuss their stellar population in the following. 

The fractional contributions ( i.e. Xj in Equation 1) of these 
stellar populations to the total stellar emission are shown in 
Figure |5] We found that they are very different for different 
types of ULIRGs. The young stellar population is more dom- 
inant for the Hll-like ULIRGs than the composite and Seyfert 
2 ULIRGs. The stellar populations of intermediate-age and 
old stars are more dominant in the Seyfert 2 ULIRGs than 
Hll-like and composite ULIRGs. The older stellar popula- 
tions and the significant power-law component in Seyfert 2 
ULIRGs suggest that they are at a more evolved stage in the 
evolution track towards AGNs. The fraction distribution of 
stellar populations of composite ULIRGs suggests that they 
are at the transitional stage between Hll-hke and Seyfert 2 
ULIRGs. 

4. PHYSICAL PARAMETERS OF DIFFERENT TYPES OF ULIRGS 

The results of stellar population analysis of ULIRGs sug- 
gest the possible evolution sequence from Hll-like, compos- 
ite, to Seyfert 2 ULIRGs. The Hll-like ULIRGs are at an ear- 
lier evolutionary stage, the composite ULIRGs are at a transi- 
tional stage, and the Seyfert 2 ULIRGs are at a more evolved 
stage. In this section, we discuss the physical parameters of 
ULIRG galaxies, namely, the mean stellar age, stellar mass, 
and gas mass, which may be related to the possible evolution- 
ary sequence. 

4.1. The mean stellar age 

Based on the results of the stellar population anal- 
ysis, we can get the l ight-weighted mean stellar age 
(icid Fernandes et al.ll2005h : 

N 

(logt,)L ^^Xj \ogtj, (2) 

1=1 



which is an indicator of star formation history. The other such 
an indicator is D„(4000) obtained by the MPA/JHU group^, 
which is defined as the average flux density ratio for the bands 
3850 A - 3950 A and 4000 A - 4100 A (iBruzual A.II1981 
iBalo gh et alJll999h . We compared them (Figure |6ll and found 
that they are well correlated, with the Spearman Rank-Order 
Correlation Coefficient of rg — 0.90. The best fitting is 

i:i„(4000) = -(0.24±0.06) + (0.18±0.01)(logi*)L. (3) 

The distributions of {logt^,)L for Hll-like, composite and 
Seyfert 2 ULIRGs are shown in the right panels of Figure |6] 
The mean values of (logt^,) l for three types of ULIRGs, HII- 
like, composite, and Seyfert 2, are 7.6, 8.0 and 8.4, with the 
standard deviations of 0.4, 0.5 and 0.4, respectively. The 
increase of these mean values is consistent with the evolu- 
tion sequence from Hll-like ULIRGs and composite ULIRGs 
to Seyfert 2 ULIRGs. The young stellar population (tj < 
10* yr) are always present in ULIRGs (see Figure|6) due to the 
enhanced star formation in galaxy merger process. The large 
deviations of {logt^)^ for these types of ULIRGs are proba- 
bly caused by the complex dynamical process and star forma- 
tion history of merger systems. Therefore, both {logt^,)L and 
D„(4000) are only the rough age indicator of ULIRGs. 

4.2. The stellar mass 

ULIRGs experience the continuing conversion of gas to 
stars. The stellar masses in ULIRGs are mostly contributed 
by intermediate-age and old stars (see /ij distribution in Fig- 
ure|4]i, and thus are a possible indicator for the ULIRG evolu- 
tion. 

The stellar mass of a ULIRG in the central 3" for the 
SDSS fiber spectrograph, A4 fiber, in unit of Mq, can be esti- 
mated from the the stellar mass parameter Aicor.tot given by 
STARLIGHT^: 

Mfiber = 7Ucor.totxl0"^^x47rd2x(3.826xl0^^)-\ (4) 

where d is the luminosity distance in cm. The aperture effect 
must be corrected b y using the SDSS Petrosian magnitude 
dBlanton et al.ll2001h . ™pctro, and fiber magnitude, mfiber, so 

' http://www.mpa-garching.mpg.de/SDSS/DR7/ 

* http://www.starlight.ufsc.br/papers/Manual_StCv04.pdf 
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Fig. 7. — Aperture-corrected stellar masses, of ULIRGs calculated 

from the STARLIGHT pai'ameters are compared with the stellar masses esti- 
mated from the SDSS photometric data by the MPA/JHU group (webpage in 
footnote No. 5). The soUd line is the best fit for all ULIRGs. The total dynami- 
cal masses of torn ULIRGs, F09039+0503, F13428-H5608, F15250-H3609and 
F15327-I-2340, are estimated from the central velocity dispersion and plotted 
as squares with error-bar The distributions of the stellar masses for each type 
of ULIRGs are shown in the lower panels. 



that the estimated mass can be a constant fraction of the to- 
tal mass of galaxies, independent of their positions and red- 
shift. The aperture correction f actor for each SDSS photo- 
metric band is roughly given by (iHopkins et al.ll2003b : 



A = 10 



-0.4(mpctro-nitibci) 



(5) 



The spectra of ULIRGs fitted for the stellar population anal- 
ysis roughly cover the u, g, r, i and z bands. We take the 
average of the correction factors of these bands weighted by 
their uncertainty Sa, 



A 



j—u,g,r,-i,z -^j j—u,g,r,t,z 



1 



(6) 



to correct the aperture effect for each ULIRG. The corrected 
stellar mass, A^*, in unit of A/0, which is a constant fraction 
of the total stellar mass in a ULIRG, then can be estimated by 



l0g(Al*) = log{M fiber) + log(^). 



(7) 



We obtained M^, for the ULIRGs in our sample. To show 
the uncertainty of such corrected mass estimates from the 
fiber spectrum of the central part of ULIRGs, we compared 
them with the total stellar masses of whole galaxies estimated 
by the MPA/JHU group from the SDSS photometric data of 
u, q, r, i and z ba nds (see the footnote No. 5 on their webpage, 
ISalim et alj2007l) . We found that they are well correlated (see 
Figure |7|l with the Spearman Rank-Order Correlation Coeffi- 
cient rg = 0.80, which means that the aperture-corrected stel- 
lar masses calculated from the STARLIGHT parameters are 
statistically fine. The data scatter in Figure |7] roughly indi- 
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Fig. 8. — The estimated gas mass from the CO measurements (the upper 
panel) and the gas fraction of total mass (gas plus stars) in 1 ULIRGs (the 
lower panel) are plotted against the aperture-corrected stellar mass M » . The 
dashed lines are the best fittings to the 9 ULIRGs except for the outlier, an 
Hll-like ULIRG F08572+3915. We also plotted the gas masses of another 
9 Hll-like ULIRGs which have CO observations available in literature but 
no SDSS spectra for stellar population analysis and hence no stellar mass. 
The mean value for Hll-like ULIRGs ( log A^, = 10.70, see Figure|7] is 
adopted for these 9 Hll-like ULIRGs. 



cates the uncertainty of the mass estimates. 

The most important feature in Figure |7] is the obviously 
different distributions of log for Hll-like, composite and 
Seyfert 2 ULIRGs, whichever the masses estimated from ei- 
ther spectral or photometric analysis are used. The mean of 
the log A^, distribution for the mass from the spectral anal- 
ysis for three types of ULIRGs, Hll-like, composite, and 
Seyfert 2, are 10.70, 11.17 and 1 1.40, with the standard devi- 
ations of 0.45, 0.38 and 0.25, respectively. Seyfert 2 ULIRGs 
always have * > 10^°^ A/q, larger than the masses in HII- 
like or composite ULIRGs (see Figure|7]i. The increase of the 
stellar masses along the sequence from Hll-like ULIRGs to 
the composite ULIRGs and Seyfert 2 ULIRGs is consistent 
with the standard evolutionary scenario of ULIRGs. 

4.3. Gas in ULIRGs 

We noticed that Seyfert 2 ULIRGs mostly have M* > 
101"-^Mq (see Figure EJ. However, many Hll-like ULIRGs 
or composite ULIRGs have much smaller stellar mass, A^^,, < 
10^'^-^ Mq. If the evolutionary sequence from Hll-like 
ULIRGs to composite ULIRGs and Seyfert ULIRGs is true, 
there should be enough gas to be converted to stars. The gas 
mass in ULIRGs is a key parameter to verify the standard evo- 
lution scenario of ULIRGs. 

We collected the CO observations of ULIRGs to study their 
gas content. Up to now, ^ 60 ULIRGs and L I RGs have mea- 
surements of CO(l-O) dSanders et all 119911: ISolomon et"!!] 
119971; IChunget all 120091) . However, onlv 10 of them have 



^ In our sample, the total dynamical masses of four ULIRGs, 
F09039-I-0503, F13428-I-5608, F1525b+3609 and F15327-I-2 340, can be cal- 
culated from the central velocity dispersion iPasvra et al.ll200 6). As sug- 
gested by the referee, we got logMdynlMQ] of 10.9, 11.1, 10.9, and 11.2 
and plotted them in the Figure|7]for comparison with the total stellar masses 
derived by the stellar population analysis. Considering the typical gas content 
of ULIRGs of ~ 10^ Af0 (see Section 4.3), we think different kinds of mass 
estimates are consistent with each other 
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the optical spectra available from SDSS, and we have got the 
relevant ma ss estimates. CO obse rvations of 9 ULIRGs were 
obsei-ved bv lSolomon et alj (Il997h usin g the IRAM 30m tele- 
scope, and 1 ULIRG by IChung et af] (2009) using FCRAO 
14m telescope. The gas mass (in unit of of AIq) can be es- 
timated from the CO lumi nosity (in unit of K km s^^ pc^) 
dPownes & Solomonll99l) . 

Mh, - 0.8 X Leo- (8) 

After the CO luminosity is converted to the same cosmol- 
ogy adopted in this paper, we got the gas masses for these 10 
ULIRGs (see Figure|8]l. They are between lO^-^ - IO^^ ^Mq, 
typical for normal gas-rich galaxies, except for the only one 
Hll-like ULIR G, F08572-H3915, which have a much sm aller 
CO luminosity ( ISolomon et 011119971: IChung et al.ll2009D and 
hence the smaller gas mass. Clearly, more gas exists in 
ULIRGs with a sm aller stellar mass. We further collected 
CO measurements (iSanders et al.lll9911:ISolomon et al.|[T997l: 
Chung et al.1 l2009l) of 9 Hll-like ULIRGs dVeilleux et al] 
19917 which have no SDSS spectra available for stellar pop- 
ulation analysis above. We take the mean log A^, for these 
Hll-Uke ULIRGs, and plotted them in Figure[|l 

The gas mass from the CO measurements (the upper panel 
of Figure|8]l and the gas fraction of total mass (the lower panel 
of Figure [S) seem to be anti-correlated with the stellar mass 
for the ULIRGs (r^ = -0.60 and = -0.88), if we do 
not consider the outlier F08572H-3915. The anti-correlation 
is only marginally significant because of the small sample of 
data available. The possible anti-correlation is preserved for 
the gas fraction and the stellar mass estimated from photomet- 
ric analysis. 

Hll-like ULIRGs with a gas fraction > 10% follow the anti- 
correlation. Gas in ULIRGs may be eventually converted to 
stars. We noticed that about half of Hll-like ULIRGs have 
a stellar mass less than 1O^°'^M0 (see Figure IT). Hll-like 
ULIRGs in Figure |8] have a gas mass less than ~ 10^^ Mq 
and a gas fraction less than ^ 15%. For those less massive 
Hll-like ULIRGs, there is not enough gas to form stars, and 
they are not massive enough in total to evolve to Seyfert 2 
ULIRGs (M^ ~ IO^A/q). Only massive Hll-like ULIRGs 
(A^* > IO^^'^Mq) may follow the standard evolutionary sce- 
nario of ULIRGs and have a possible connection between the 
starbursts and AGNs. 

5. DISCUSSIONS 

Stellar population analysis for narrow-line ULIRGs show 
the systematical changes of stellar age and mass along the se- 
quence from the Hll-like ULIRGs, to the composite ULIRGs 
and Seyfert 2 ULIRGs. Here we discuss several factors which 
may affect the results of stellar population analysis. 

5.1. Power-law AGN contribution for spectral continuum 

During the stellar population analysis with STARLIGHT, 
we used a power-law component, F{\) oc A", with a fixed 
a = — 1 . 5 to represent the AGN contribution to the continuum 
spectrum (ICid Femandes et al. 2004b, 2005). The AGN con- 
tribution is significant for composite, Seyfert 2 and LINER 
ULIRGs, but not for Hll-like ULIGRs. However, the index 
for the power-law contribution, a, could be —2.0 to —1.0 (see 
iNatah et"aLlll998 ). To evaluate the influence of the power-law 
index a in the stellar population analysis, we re-do the fittings 
to the all spectra of 160 ULIRG sample with different values 
of a, from —1.0, —1.25, —1.5, —1.75, to —2.0, for ten times 
each with different random seeds. We found that the change 
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Fig. 9. — Same as Figure|6]but only for ULIRGs in the redshift range of 
0.20-0.24. 



of by different a values is less than ^1%, which is not ad- 
equate to use to judge the best a for fitting. Changes of 
the stellar mass and the mean stellar age caused by different a 
are always less than 10%. Therefore, the adopted power-law 
index of a — —1.5 is fine for the stellar population analysis. 

5.2. Different bases of simple stellar populations 

The evolutionary population synthesis method is mainly 
based on the stellar evolution model. We have used a base 
of 45 simple stellar populations (SSPs) for the stellar popu- 
lation analysis. The different bases of SSPs may affect the 
fitting results. Following the same procedures described in 
Section 3, we re-do the fittings to the spectra for 160 ULIRGs 
with a base of 150 SSPs (25 ages and 6 metallicities) and a 
base of 15 SSPs ( 15 ages and 1 metallicity) extracted from 
iBruzual & Chario t! ^003). Although the fractional contribu- 
tion of each stellar population component does not remain the 
same and may vary about 20%, the mean stellar age and stel- 
lar mass do not change significantly, less than ^10%. Their 
distributions are consistent with those shown in Section 4 for 
different types of ULIRGs. 

5.3. The aperture effect 

We made stellar population analysis for the SDSS spectrum 
of the central 3" region. The aperture effect has been cor- 
rected for the mass estimates, but not for the mean stellar age. 
To evaluate the aperture effect, we divided the ULIRG sam- 
ple to several redshift bins, and re-do the analysis. In each 
redshift bin, the linear size of the aperture is almost the same 
for all ULIRGs, and the aperture effect is almost the same. 
Here we take the redshift bin of 0.20 < z < 0.24 as the ex- 
ample which has the largest number of ULIRGs in the bin. 
As shown in Figure |9] the correlation between £'„(4000) and 
{logtf) L keeps the same as the whole sample of ULIRGs. The 
distributions of {logt^^L for Hll-like, composite and Seyfert 
2 ULIRGs in the right panel of Figure |9] also show similar 
offsets as the whole sample. So do for other redshift bins. 
Therefore, we believe that the aperture effect does not influ- 
ence our results. 

5.4. Extinction 

In the stellar population synthesis model of STARLIGHT 
(ICid Femandes et al.l 2005: Mateus et al. 2006), the extinction 
was considered as a foreground screen and fitted by one pa- 
rameter, the y-band extinction Ay- But the young, interme- 
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Fig. 10. — Distribution of tlie extinction pai'ameter Ay given by tlie 
STARLIGHT model for Hll-like ULIRGs, composite ULIRGs and Seyfert 
2 ULIRGs. 



diate and old stellar populations have different extinctions in a 
galaxy. Particularly the young stellar population was more ex- 
tinguished than the others, being in places very dust obscured. 
Different Ay should be assumed for the different populations. 
However, this can not be achieved in the model. 

From the stellar population analysis to ULIRGs, the 
distributions of extinction parameter Ay given by the 
STARLIGHT model for Hll-like ULIRGs, composite 
ULIRGs and Seyfert 2 ULIRGs are very similar, as shown 
in Figure [10] The means of Ay for three types of ULIRGs 
are 0.91, L09 and 0.92 with the standard deviations of 0.28, 
0.39 and 0.36, respectively. The extinction parameters Ay 
of Hll-like and Seyfert 2 ULIRGs on average are slightly 
smal ler than that of comp osite ULIRGs, which is consistent 
with lVeilleux et al.1 ( l2009l) . 

The distributions of the mean stellar age for different 
ULIRGs discussed in Section 4 are probably not influenced by 
the simplified treatment of extinction in the model. The mean 
stellar age given by the stellar population synthesis model are 
well correlated with results of D„(4000) (see Figure|6]l, which 
verifies the mean age. On the other hand, in ULIRGs, the stel- 
lar mass is dominant by old stellar populations (see the right 
panels of Figure lU, which are expected to suffer less from 
extinction than young stellar populations. Therefore, the sim- 
plified extinction in the model probably does not affect our 
conclusions for ULIRGs. 

6. CONCLUSIONS 

We analyse the stellar populations for a sample of 160 
narrow-line ULIRGs. They are optically classified as 32 HII- 
like ULIRGs, 74 Seyfert-HII composite ULIRGs, 6 LINER 



ULIRGs and 48 Seyfert 2 ULIRGs. We found fliat that along 
the sequence of Hll-like, composite and Seyfert 2 ULIRGs, 
both the mean stellar age and the aperture-corrected stellar 
mass increase. This supports the standard evolutionary sce- 
nario of ULIRGs in which the Seyfert 2 ULIRGs are in a late 
stage of ULIRG evolution with a large mean stellar age and 
a large steflai- mass (-- IG^A/q) while Hll-like ULIRGs are 
in an early stage with a young stellar age and a small stellar 
mass on average. 

Whether Hll-like ULIRGs have enough gas for starburst, 
so that they can evolve to Seyfert 2 ULIRGs? We collected 
CO measurements of 10 ULIRGs which have the SDSS spec- 
tra for the stellar population analysis. The gas mass fractions 
seem to be anti-correlated with the stellar masses for massive 
ULIRGs. All ULIRGs in our sample with CO measurements 
have a gas mass less than ^ Hll-like ULIRGs with a 

small stellar mass (A^* < IO^^'^Mq) do not possess enough 
gases for starburst, and therefore have no evolution connec- 
tions with massive Seyfert 2 ULIRGs. We conclude that only 
massive Hll-like ULIRGs and composite ULIRGs may follow 
the evolution sequence toward AGNs. 
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